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Abstract 

We study event-by-event dynamical fluctuations of various particle ratios at different ener- 
gies. We assume that particle production in final state is due to chemical equilibrium processes. 
We compare results from resonance gas model with available experimental data. At SPS ener- 
gies, the model can very well reproduce the experimentally measured fluctuations. We make 
predictions for dynamical fluctuations of strangeness and non-strangeness particle ratios. We 
found that the energy-dependence is non-monotonic. Furthermore, we found that fluctuations 
strongly depend on particle ratios. 

1 Introduction 

Understanding dynamical properties of hot and dense matter is a key question in heavy-ion coUi- 
sions experiments. Event-by-event fluctuations have been suggested ^ElE] to provide comprehen- 
sive characteristics of this matter. They are crucial observations to check hypothesis of chemical 
equilibrium in heavy- ion collisions jl] . Event-by-event fluctuations of very few particle ratios have 
been studied in several experiments at SPS and RHIC energies [SIIHIE]. In this letter, we pur- 
pose to study energy-dependence of event-by-event fluctuations in hadron resonance gas model. 
Hypothetical phase transition to and from quark-gluon plasma can be characterized by large fluc- 
tuations in particle yields ^ElE], which are accompanied by volume fluctuations. The latter 
can be eliminated, when particle ratios are considered In this letter, we try to give answers 
to the questions, whether strangeness quarks enhance dynamical fluctuations and whether critical 
endpoint can be localized by means of dynamical fluctuations. Since experimental measurements 
of event-by-event fluctuations are very much limited, we hope that our predictions encourage spec- 
ifying certain particle ratios and measure their dynamical fluctuations. In this letter, we make 
predictions for dynamical and statistical fluctuations of different particle ratios in dependence on 
energy. The values of fluctuations strongly depend on particle ratios. In some particle ratios, 
dynamical fluctuations are smaller than statistical ones. In other particle ratios, dynamical fluc- 
tuations are slightly greater than statistical ones. We also found that fluctuations of hybrid and 
cascade baryons are dominant. All these predictions are phenomenologically of great interest. 



2 Model 

Pressure in hadronic phase is given by contributions from all hadron resonances treated as a free 
gas [HIIHIEIE]- The resulting pressure accounts for free as well as for strong interactions between 
resonances. In previous work |llj. we have proved that the thermodynamics of strongly interacting 
system of hadron resonances can be approximated by an ideal gas of stable and resonance hadrons. 

At flnite temperature T, strangeness /is and iso-spin and baryo-chemical potential /ts, 
pressure of one hadron reads 

1± 7ABAsA/3e^ , (1) 
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where s{k) = (fc^ + to^)^/^ is single-particle energy and ± stands for bosons and fcrmions, respec- 
tively, g is spin-isospin degeneracy factor. 7 = 7^7^* are quark phase space occupancy parameters, 
where n and m being number of light and strange quarks, respectively. In this letter, we explicitly 
use the equilibrium value for 7, the unity. A = exp(/i/T) is fugacity, where /i is chemical potential 
multiplied by corresponding charge. 

Quark chemistry is given by relating hadronic chemical potentials to quark constituents, fis — 
and fis = fJ-q — fJ-s, where q and s being light and strange quark quantum number, respectively. 
The baryo-chemical potential for light quarks is /Zq = {fiu + fj.d)/'^- fJ-s is calculated as a function 
of T and fis under the condition of strangeness conservation. Iso-spin chemical potential fii^ = 
(m« - A*<i)/2. 

Particle number density is given by derivative of partition function in Eq.^with respect to the 
chemical potential of interest. Fluctuations in particle number are given by susceptibility density, 
which is second derivative with respect to chemical potential. 

After freeze out, hadron resonances decay either to stable particles or to other resonances. 
Particle number and fluctuation density in final state have to take into account this chemical 
process. 

(nf"'^') = (nf— *)+^6,_K), (4) 

((An,_0') = 6,^,;(l-6,-..)S>+^^,%((An,r) (5) 

where bj^i being branching ratio for decay of j-th resonance to «-th particle. Chemical freeze out 
is characterized by s/T^ where s is the entropy density. 
Fluctuations of particle ratio 71-1/712 are [3] 



((Ani)^) {(An,)^) {An, An,) 



which include dynamical as well as statistical fluctuations. Third term of Eq.El counts for fluctu- 
ations from hadron resonances that decay into particle 1 and particle 2, simultaneously. In such 
a mixing channel, all correlations including quantum statistics ones are taken into account. Obvi- 
ously, this decay channel results in strong correlated particles. To extract statistical fluctuation, 
we apply Poisson scaling in mixed decay channels ^, 

' (ni) {n2) 

Subtracting Eq. 13 from Eq. we get dynamical fluctuations of particle ratio ni/n,. 

{"ni/noJdyn ~ / \2'/\2 /\/\ ^ ' 



3 Results 



For the first time, experimental measurements of dynamical fluctuations of particle ratios are 
systematically confronted with theoretical predictions in Fig. ^ There was an earlier attempt 
to compare with preliminary results for K/tt reported in 3 . Is has been found that theoretical 
and experimental ratios of dynamical to statistical fluctuations are compatible with each other. 
Individual fluctuations themselves were not. 

^Experimentally, there are various methods to construct statistical fluctuations Frequently used method is 
the one that measures particle ratios from mixing events. 
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Hadron resonance gas model is obviously able to reproduce all measurements. For K'^ /t:^ ratio, 
experimental data covers a wide range of square root of center of mass energy ^/s . While for 
(p+p)/(7r^ +7r~) ratio, experimental data available so far has been measured at SPS energies only. 
That the results from hadron resonance gas model agree well with the existing experimental data, 
allows us to make predictions for other particle ratios as will be given in Fig. |21 On other hand, 
we can now use the model to systematically investigate dependence of event-by-event fluctuations 
on ^/s. Few comments are in order at this moment. 

• Dependence of event-by-event fluctuations of particle ratios on ^/s is non-monotonic 

• Fluctuations can be suppressed and/or enhanced at different ^/s 

• Strangeness fluctuations are positive and enhanced with ^/s. Anti-strangeness fluctuations 
are smaller than strangeness ones. There are remarkable minima at top SPS energies 

• Non-strangeness fluctuations are negative, i.e., statistical fluctuations are greater than dy- 
namical ones, especially at low energies. 

• In general, event-by-event fluctuations at high energy increase with increasing energy. There 
is only one exception from this empirical role; K~ /■k~ ratio. Fluctuations of K~ /■k~ ratio 
exponentially decrease with ^/s. 

In left panel of Fig.^ we compare experimentally measured fluctuations of /tt^ ratios with 
resonance gas model results. There is a good agreement at SPS energies. At RHIC, measured fluc- 
tuations are slightly above the theoretical line. The explanation for this disagreement is twofolds. 
First, RHIC measurements are still preliminary |2j. Second, we refer to our previous study of 
particle ratios in heavy ion collisions ^Hl- We have concluded that thermal models with 7 = 1 
slightly overestimate particle ratios at RHIC. To reproduce particle ratios, precisely, one has to 
allow 7 to take values different from unity Inserting large (n), averaged particle number density, 
in Eq. |S1 apparently results in small a. 

In right panel in Fig. ^ dynamical fluctuations of [p + p)/(7r+ -I- 7r~) ratios are depicted as a 
function of ^/s. According to Eq. |H1 negative values of dynamical fluctuations are to be understood 
as a reason of much dominant statistical fluctuations. Nevertheless, using our model, we now have 
a quantitative estimation of both dynamical and statistical fluctuations. At low SPS energy, there 
is a good agreement between our model and the experimental results. 

A much more important quantity is the ratio of dynamical fluctuations to statistical ones. This 
quantity will be given in Tab.Qlat specific energies. Further advantage of right panel of Fig.^is to 
illustrate the ability of our model. To author's knowledge, there is no another experimental data 
available. 



Our predictions for dynamical fluctuations of different particle ratios are depicted in Fig. |21 
Fluctuations of strangeness particle ratios are given in left panel. Right panel shows fluctuations of 
non-strangeness particle ratios. Depending on particle ratios, some fluctuations decrease and others 
increase with y/s. As we discussed above, the energy dependence is apparently non-monotonic. 
Some particle ratios have negative dynamical fluctuations. Fluctuations of other particle ratios 
move from negative to positive values. 

Negative fluctuations are explicitly obtained from non-strangeness particle ratios; dynamical 
fluctuations of ^/tt", p/7r+ and 7r+/7r~ ratios are overall negative. Statistical fluctuations of such 
particle ratios are much greater than dynamical fluctuations. The fluctuations of p/tt~ and njix^ 
have a remarkable dependence on At low energy, i.e., AGS and SPS energies, they are negative. 
Within a relatively short energy interval, their values switch to positive sign. At higher energies, 
they monotonically raise with increasing -y/s. 

The values of dynamical fluctuations of strangeness particle ratios are greater than that of non- 
strangeness ones. Their dependence on -^s is also non-monotonic. While fluctuations from /'k~ 
exponentially decrease with increasing y^, fluctuations from A/tt and S^/vr have minimum values 
at ~ 10 GeV. 

■^Since available data is still preliminary, we can, for this moment, restrict our calculations to equilibrium value. 
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Fig. 1: Left panel: dynamical fluctuations of /v+/7r+ ratio as a function of square root of center of 
mass energy (curve) compared with experimental results (circles). Results from SPS experiments ^ 
are drawn in open symbols. RHIC results are still preliminary (solid circles). Right panel shows 
dynamical fluctuations of (p + p)/(7r+ +7r~) ratios. Negative values are an indication to dynamical 
fluctuations that are smaller than statistical ones. 
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Fig. 2: Predictions for dynamical fluctuations of strangeness (left panel) and non-strangeness (right 
panel) particle ratios. Values of dynamical fluctuations at different energies strongly depend on 
particle ratios. 



Comparing strangeness with non-strangeness fluctuations, we find that replacing pion by its 
anti-particle has almost no influence on strangeness dynamical fluctuations. Non-strangeness dy- 
namical fluctuations arc dramatically changed, when pion in denominator ^ has been replaced by 
its anti-particle. For instance, for p/i: ratio 

((Anp/,+ )') ^ ((An„)^)+2((An,)2) (9) 
((Anp/,-)') ^ 3((An„)^) (10) 

As in right panel in Fig. El while fluctuations of first ratio move from negative to positive values, 
second ratio remains negative at all energies. In hadronic phase, i.e. particle ratios, quarks are 
strongly confined into hadronic states. It is believed that quarks in the unconfined phase may be 
strongly correlated jj^. It would be interesting to verify above expressions. In doing this, we have 
to take into consideration volume fluctuations on lattice. Resonance gas model can not be applied 
at temperatures higher than critical one. 

^Assuming that particle ratios can mathematically be seen as fractions. 
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4 Discussions and conclusions 



We find that the values of dynamical fluctuations depend on particle yields. According to QCD, 
hadronic matter under extreme conditions of high temperature and pressure is conjectured to have 
dynamical transitions into partonic matter. Assuming that we put certain particle ratio under these 
conditions, we want to see the change in its dynamical fluctuation before and after transition. This 
is equivalent to our numerical study at different energies. We start with hadronic matter of two 
particles and end with partonic mater, Eq. Inland Eg. 1101 

Fluctuations of quark number have been studied in lattice QCD [El ■ It has been found that 
temperature-dependence (equivalent to ^/s-dependence) of fluctuations is dominated by analytic 
part of partition function. Across critical temperature, there is a smooth increase of fluctuations 
with increasing temperature. For our analysis, critical temperature is not exactly specified. Our 
results on fluctuations of individual particle are not shown here. Nevertheless, we can make 
following statements: We find that fluctuations of mesons (individual particles) have almost same 
structure as lattice fluctuations and that baryons have maximum value at -^/s ~ 10 GeV. At higher 
energies, baryon fluctuations decrease. Quantitative comparison with lattice is - of course - not 
possible. As discussed in previous section, thermodynamics in hadronic phase can be reproduced 
when taking into account many stable and resonance hadrons in the partition function. A free gas 
of one hadron is not able to reproduce the quantitative thermodynamics below critical temperature. 

In Tab^ values of <Jdyn/'^stat quantity of different particle ratios are given at specific energies 
related to certain beam energies at top SPS and RHIC. It would - of course - be of great interest 
to verify these values experimentally. Our predictions specify the values of dynamical fluctuations 
related to statistical ones in all heavy ion collisions experiments. Experimentalists can now decide 
which particle ratios shall be measured. It depends - among others - on how large are dynamical 
fluctuations compared to statistical ones. 
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Tab. 1: Ratios of crdyn/o'stat from different particle ratios at specific energies related to known beam 
energies at SPS and RHIC. With this quantity, we estimate how large are dynamical fluctuations 
compared to statistical ones at chemical freeze out in different heavy ion collisions experiments. 
Validity of these values depend on validity of chemical equilibrium freeze out scenario. 

Strangeness fluctuation are given in right panels of Fig. ^ and Fig. [3 The dependence of 
strangeness fluctuations on energy is also non-monotonic. For /tt^ ratios, fluctuations have 
a negative dependence on energy at SPS. There is a minimum around ^/s ~ 20 GeV. At higher 
energies, fluctuations smoothly increase. Fluctuations of A/tt and S'^/tt ratios have also mini- 
mum values at ^/s ~ 10 GcV. K~ /■k~ ratio has a completely different behavior. Its dynamical 
fluctuations of this particle ratio exponentially decrease with energy. 

From our predictions for non-strangeness fluctuations drawn in right panel in Fig. |21 we find 
that event-by-event dynamical fiuctuations depend on particle ratios. There is a rapid increase 
at low energy. At high energies, fluctuations smoothly increase. There is a remarkable sharp in- 
crease in p/tt" and n/7r+ ratios within a short range of energies. The values of energies are different. 

In final conclusion, we propose to study event-by-event dynamical fluctuations of different par- 
ticle ratios in heavy ion collisions experiments. We have shown that baryon to meson fluctuations 
are much larger than meson to meson ones. Fluctuations of hybrid and cascade baryons are much 
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larger than the statistical ones. According to our results, we now have a framework to study 
fluctuations of particle ratios, systematically. As mentioned above, fluctuations of particle ratios 
eliminates volume fluctuations. The latter are still included in lattice simulations. 

We assumed that particle production is due to chemical equilibrium processes in final state, 
i.e., 7 = 1. That our models can very well reproduce experimental measurements means that the 
equilibrium freeze out scenario is apparently proved, especially at SPS energies. At RHIC, precise 
measurements are needed. Energy scan down to ^/s = 10 GeV turns to be a crucial step to verify 
the worthwhile behavior of particle production jl2lll8j and now dynamical fluctuations. 

According to this systematic study, it was not possible to point out certain energy or region 
of energy, at which dynamical fluctuations sharply increase. The speculations on manifestation 
of critical endpoint by a rapidly increase of dynamical fluctuations of particle ratios can not be 
verified by this model. Furthermore, the region of unconfined phase can not be specified, precisely. 

It would be interesting to extract information about the role of different decay channel in the 
energy-dependence of dynamical fluctuations. It will be a further propose to study the effect of 
chemical non-equilibrium processes of event-by-event dynamical fluctuations of particle ratios. 
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